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Abstract 

Sitnikov problem, consisting two close binaries and a third small body 
is considered, leading to a rapid ejection of the small body from the bi- 
naries. This mechanism is proposed as an explanation of jets in many 
astrophysical systems. Choosing appropriate initial condition relativistic 
final velocities can be achieved. 

1 Introduction 

Sitnikov problem [T] is a special 3-body problem in which two massive 
bodies follow keplerian orbits and a third, with infinitesmal mass, travel 
along straight line perpendicular to the primaries' plane, through their 
barycentre. This problem acquired considerable attention of physicists 
mostly in the case of oscillatory motion [2]- Here, however, we focus on 
unbound trajectories possible in this system and propose it as a source of 
jets observable in many stellar and galactic systems. 

One of the most striking observable effect in the Universe is a produc- 
tion of jets in many astrophysical systems which can stretch on a large 
distances (even million of light years) and can posses relativistic speeds. 
The usual model for production of these jets is a single rotating black 
hole surronded by accretion disk with strong magnetic field [3|. However, 
also differnet models are considered [4]. Here, we propose a following 
mechanism: two large masses (possibly black holes) forming closed binary 
systems, with high eccentricity elliptical orbits, surronded by accretion 
disk but with no magnetic field needed. Black holes can be of stellar 
masses as well as supermassive in the center of the galaxies. 

Let's consider a small particle at a distance zo from the binaries, per- 
pendicular to their plane. Assume that the binaries have the same mass 
and that they orbits are ellipses with high-eccentricity. Neglect any rel- 
ativistic effects and treat the problem purely Newtonian. Because of the 
gravitational attraction from the primaries, the particle will follow straight 
line along the z direction. We have got here the special case of the fa- 
mous restricted three-body problem which is non-integrable system. If 
at the time of particle crossing the center of mass of the binaries they 
are in their periastron position (denoted B on the Fig. 1) nothing inter- 
esting happens - the particle will just follow usually oscillatory motion. 
To exclude this possibility let's assume that the binaries are then not in 
their periastron position. In this case the situation is more interesting 
and choosing appropiately the position of the primaries in the moment of 



maximum encounter, energy can be transferred from them to the particle. 
We are interested in the question: what is the maximum ejection veloc- 
ity of the particle from the closed binary-system? This problem is related 
to the gravitational assist, commonly used in Solar System exploration [5]. 

We expect the maximum effect in the case where the binaries are 
coming to their periastron position just before the particle attain their 
center of mass. In this case, when all three bodies are minimally separated 
from each other and the particle is still before the center of mass point, it 
will get a big push forward ejecting it from the system, which will not be 
cancelled by later slowing down of the particle. 



2 Numerical computation 

To answer the question we performed numerical integration choosing the 
initial position of the binaries in their periastron position (denoted B 
on the diagram), the initial velocity of the particle as zero and we were 
changing the initial distance zo from the particle to the center of the mass 
of the system as to find the best initial position leading to a maximal 
ejection velocity of the particle along a straight line away from the system. 
The two masses forming binary system follow keplerian orbits and to 
determine their position in any time we used the Kepler equation. Our 
system of equation is: 
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where u is calculated at any time t from the Kepler equation [6l [7] : 

= u — e ■ sin{u) (2) 



The trajectory of the particle is integrated using the standard 4th or- 
der Runge-Kutta method. The code we used is presented in the appendix. 

We assumed at first the following data of the closed binary system: 
semimajor axis a = 6 • 10 7 m, semiminor axis b — 2 • 10 7 m, G (gravita- 
tional constant) and M (mass of the star) product G ■ M — 10 22 . We've 
found the maximum velocity which the particle received as 1.31 ■ 10 s — 
when crossing the plane of the binary system. Later the gravitational 
attraction was slowing down the particle but assymptotically (far away 
from the system) it still kept the velocity equal to 0.73- 10 7 ™ that is above 
20% the speed of light. This happened when the initial distance zq was 
equal to —4.931 • a. 

As a second example we considered the PSR1913+16-like system and 
choosed the following data: semimajor axis a = 1.95 • 10 9 m, eccentricity 
e = 0.617131, orbital period = 27900s, GM product G • M = 2 • 10 20 . 
In this case because of the smaller masses and larger parameters of the 
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Figure 1: Sitnikov problem, two large masses following keplerian orbits and a 
small particle which can attain large velocity away from the system. 



ellipses the effect was much smaller but still the particle attained 538 — 
final velocity for the initial distance zo = —3.25 • a. 



3 Conclusion 

It is almost 350 years since Newton formulated his famous laws of motion 
and gravity. He quickly showed how to solve they for two-body problem 
but his attempts to solve three-body problem failed. Since this time the 
brightest men tried to solve it but they were unsuccesful too. Apart from 
some special solution (like Lagrange one, 1772) no general exist [6] [7J [8] . 
In this paper we proposed a special 3-body case known as a Sitnikov 
problem as an explanation of jets in astrophysical systems. This effect 
can also be used in future interstellar travels, Sitnikov problem playing a 
role of a gravitational assist. 



4 Appendix 

The program was written using Quincy2005 (MinGW gcc 4.2.1) and exe- 
cuted under Windows XP on Pentium M processor (IBM ThinkPad X40 
machine). To save the results in a file run: program. exe>sitnikov.txt 

#include<iostream> 

#in elude <cmath> 

using namespace std ; 

const double pi =3. 1415926; 

double x , r , z , v , a , b ,GM, okres , e , t , dt , 

klz , k2z , k3z , k4z , klv , k2v , k3v , k4v ; 

double keplereq (double t); 

int main() 
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{a=6e + 07; b=2e + 07; 
GM=le + 22; //each mass 
e=sqrt (1 — b*b/(a*a) ) ; // eccentricity 
okr es=2* pi *pow (a, 1.5) / sqrt (2 *GM) ; 
dt=0.01; t=0; z = -4.931*a; v = 0; r=a*(l-e)/2; 

cout «"a=„"«a«" „jGM=„"«GEMC<" „„e=„"«e«" „„zO=„"«z«"\n\n" ; 

for (t=0; t< = 100; t=t+dt ) 
{ klz=v* dt ; 
klv=— 2<3VI* z /pow ( ( z* z+r *r) , 1 . 5 ) * dt ; 
k2z = (v + 0.5*klv)* dt ; 

k2v=-2=iGM*(z + 0.5*klz)/pow( ( (z + 0.5*klz)*(z + 0.5*klz) + r*r),1.5)*dt; 
k3z = (v + 0.5*k2v)*dt ; 

k3v=-2=fGM*(z + 0.5*k2z)/pow( ( (z + 0.5*k2z)*(z + 0.5*k2z) + r*r),1.5)*dt; 
k4z = (v + 0.5*klv)* dt ; 

k4v=-2<M* ( z+k3z )/pow( ( ( z+k3z )* ( z+k3z)+r*r ) , 1 .5)* dt ; 

z=z + (klz+2*k2z+2*k3z+k4z ) /6; 
v =v+(klv+2*k2v+2*k3v+k4v ) /6; 
r=a* (1 — c*cos(keplereq(t ))) / 2 ; 

cout « " t= J'«t«" z=„"«z«" v=„"«v«" r=„"«r«" \n" ; 

} 

} 

double keplercq (double t) 
{ double f , fmin , u , umin ; 
fmin = 100; f=0; u = 0; umin=0; 
if ( t >=o k r e s ) 

{while (t>=okres) { t=t — okres ; } } 
for (u=0;u<2*pi ; u=u + 0.01) 

{ f=u— c* sin(u) — 2*pi*t/okres ; 
if (fabs(f) < fabs (fmin)) 
{ fmin=f ; umin=u ; } 

} 

return umin; 
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